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The list of possible commensurate phases obtained from the parent tetragonal phase of 
Ruddlesden-Popper systems, A„+iB„C3„+i for general n due to a single phase transition involving 
the reorienting of octahedra of C (oxygen) ions is reexamined using a Landau expansion. This 
expansion allows for the nonhnearity of the octahedral rotations and the rotation-strain couphng. 
It is found that most structures allowed by symmetry are inconsistent with the constraint of rigid 
octahedra which dictates the form of the quartic terms in the Landau free energy. For A2BC4 our 
analysis allows only 10 (see Table III) of the 41 structures listed by Hatch et al. which are allowed 
by general symmetry arguments. The symmetry of rotations for RP systems with n > 2 is clarified. 
Our list of possible structures in Table VII excludes many structures allowed in previous studies. 

PACS numbers: 61.50.Ks,61.66.-f,63.20.-e,76.50.+g 



I. INTRODUCTION 

The Ruddlesden-Popper (RP) compounds [ij are 
layered perovskites having the chemical formula 
A„+iB„C3„+i and which consist of two slabs of BCg 
octahedra per conventional unit cell. Each slab consists 
of n layers of corner sharing octahedra of F's or O's. 
These systems either are or can be considered to be 
developed via one or more structural transitions from 
the high symmetry tetragonal parent structure shown 
in Fig. [T]for the cases exemplified by K2MgF4 {n = 1) 
and Ca3Mn2 07 (n = 2). We will refer to the RP system 
with 71 = 1 as RP214 and to that with n = 2 as RP327. 

The RP systems exhibit many interesting technological 
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FIG. 1: (Color online) The high symmetry (body centered 
tetragonal) parent structure of A2BO4 (left) and A3B2O7 
(right). The green squares are A ions. The B ions are at 
the centers of the oxygen (blue dots) octahedra. 



properties such as high Tc superconductivity [1], colos- 
sal magnetoresistance,Q metal insulator transitions, [3] 
and coupled ferroelectric and magnetic order. [^-0] Many 
of these properties depend sensitively on the distortions 
from the ideal tetragonal I4/mmm structure (see Fig. [T]) 
of space group #139 [numbering of space groups is from 
Ref. Q which appear at structural phase transitions. [Qj- 
ITsj Accordingly, the accurate characterization of their 
structure is essential to reach a detailed understanding 
of their properties. Such an understanding can poten- 
tially lead to the fabrication of new systems with en- 
hanced desired properties. It is therefore not surprising 
that one of the most celebrated theoretical problems in 
crystallography is to list the possible structures that can 
result from a single structural phase transition in which 
the (usually oxygen) octahedra are cooperatively reori- 
ented under the constraint that they are only distorted 
weakly (in a sense made precise below) . One of the earli- 
est works to address this question was that of Glazer 14 
who analyzed possible structural distortion from the cu- 
bic parent structure of CaTiOs. It turned out that a few 
of the structures he found did not actually satisfy the 
constraint of not distorting the octahedra. [Tsj For the 
RP214 systems the two principal approaches to this prob- 
lem which have been used are a) a direct enumeration of 
likely structures [l6| and b) the use symmetry. [l7| This 
last approach utilizes a very useful computer program [l^ 
to generate the isotropy subgroups of Ref. [3 In this way 
Hatch et al. fv^ gave a listing for the RP214 structure of 
possible phase transitions involving distortions at various 
high symmetry wave vectors. This listing was shown to 
be consistent with the revised results of method al.frzj 
This important work has stood unquestioned for over a 
decade. [201 H ere we show that most of structures listed 
in Refs. IITI and [l^ for the RP214 systems are incon- 
sistent with the constraint of rigid octahedra. To imple- 
ment this constraint, we assume that the spring constants 
for distortion of the octahedra are larger than the other 
spring constants of the lattice by a factor of A. Most of 
our results are obtained to leading order in 1/A, which 
we regard as an expansion parameter. This constraint 



causes the quartic terms in in the Landau free energy 
to assume a form which is less general than allowed by 
symmetry, [m In some cases this constraint causes us to 
reject structures which have undistorted sublattices, a 
situation which is counterintuitive, since it is analogous 
to having a magnetic system simultaneously having or- 
dered and disordered sublattices. Even for structures we 
our analysis allows, it is inevitable that in the structural 
phase transition the octahedra will undergo small (of or- 
der 1/A) distortions, which are observed. [2^ 

Briefly this paper is organized as follows. In Sec. II 
we enumerate the high symmetry wave vectors of the 
distortions we will consider and we discuss the role of 
the quartic terms in the Landau expansion in determin- 
ing the detailed nature of the distortions. Here we also 
develop the nonlinear constraint induced by the rigidity 
of the octahedra. In Sec. Ill we apply these ideas to 
enumerate the possible structures which are allowed via 
a single phase transition involving a distortion at these 
high symmetry wave vectors for the RP214 structure. In 
Sec. IV we extend the treatment to the analogous RP327 
(A3B2C7) bilayer system. In Sec. V we use our results 
for n = 1 and n — 2 RP systems to obtain results for 
the RP systems A„_|_iB„C3„-|-i consisting of n-layer slabs 
(with n finite) . In Sec. VI we discuss and summarize our 
results. 




FIG. 2: (Color online) The first Brillouin zone for RP sys- 
tems. There are two inequivalent X = (1/2,1/2,0) points, 
four inequivalent N points, and two inequivalent P points. 
(Wave vectors are "equivalent" if their difference involves an 
integer number of reciprocal lattice vectors G.) 



P = (1/2, 1/2, 1/2). [13 These vectors are on faces of the 
first Brillouin zone as shown in Fig. [21 The reciprocal 
lattice vectors are 

Gi = (-1,1,1), G2 = (1,-1,1), 

G3 = (1,1,-1) , (3) 



II. GENERAL PRINCIPLES 

A. OVERVIEW 

We will analyze possible distortions from the parent 
tetragonal system using a Landau-like formulation in 
which we write the free energy F as 

F - ^J2'^kAT)XkXi + OiX^) , (1) 

k,l 

where Xk is a component of an ionic displacement. A 
structural phase transition occurs at a temperature Tq 
when an eigenvalue of the matrix A becomes zero. (For 
T > Tq all the eigenvalues of A are positive.) If the zero 
(critical) eigenvalue is A^-fold degenerate, then for T near 
To one has 

N 

F^^{T~T,)Y,Ql, (2) 

/c=l 

where is the amplitude of the fcth linear combination 
of X's given by the fcth critical eigenvector of A. As 
we shall see, higher order (in Q) corrections to the free 
energy in the cases of interest involve only even powers 
of the Q's. 

As is customary (e. g. see Ref. [13), we will restrict 
attention to the cases when Q is a superposition of dis- 
placements associated with the star of the high symmetry 
wave vectors X = (1/2, 1/2,0), of N = (1/2, 0, 1/2), or of 



Instead of dealing with irreducible representations (ir- 
reps), we will develop the free energy for the most gen- 
eral structure which can be constructed using the angular 
distortions at the wave vectors of the star of X, N, or P. 

To see how the form of the fourth order potential af- 
fects possible structural distortions, consider a system 
with two order parameters Qi and Q2 related by sym- 
metry, for which the free energy assumes the form 

F = {T-To)[Ql + Ql]+u[Ql + Qlf+vQlQl{A) 

up to fourth order in Q with u > 0. As the tempera- 
ture is lowered through the value Tq the nature of the 
ordering depends on the sign of v. (See the phase dia- 
gram of Fig. [3]) If w is positive, then ordering has either 
Qi or Q2 zero. If — 4u < v < 0, ordering occurs with 
IQil — \Q2\- At the multicritical point[24,25| where v 
is zero (and also a similar sixth order anisotropy van- 
ishes) one can have ordering in an arbitrary direction of 
order parameter {Q\-Q2) space. Alternatively, analysis 
of terms of higher order than shows that in extreme 
limits ordering can occur in an arbitrary direction in Qi- 
Q2 space. Note that if we invoke only the symmetry 
properties of the system, there is no constraint on the 
parameters u and v, in which case the analysis of Ref. 
llTI would apply. However, as we will explain below, the 
picture of the lattice as consisting of oxygen octahedra 
implies a special form of the quartic terms with u > 
and V — —2m, so that in most cases only structures with 
IQil = IQ2I will actually occur. 
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FIG. 3: (Color online) The mean- field phase diagram psj 
(schematic) for the free energy of Eq. Q, showing the mul- 
ticritical point (MCP) at t; = 0. The dashed line indicates a 
continuous phase transition and the solid line a discontinuous 
one. 



B. Octahedral Constraint 

We now discuss the possible distorted configurations 
as constrained by the rigidity of the oxygen (or C ion) 
octahedra. In describing the distortions we use a no- 
tation similar to that of Ref. [13 and [13 which we deem 
more convenient for the RP systems than the widely used 
Glazer notation[l3| for pseudocubic perovskites. For a 
single RP layer, we will first consider the rotation of oc- 
tahedra through an angle 29 about the tetragonal axis. 
Giving the angle of rotation of on octahedron at the ori- 
gin fixes the rotation angles of all other octahedra in this 
layer. In the Glazer notation this would be specified as 
aPaPb'^ or aPaPb~ , where the superscripts indicate zero 
rotation about the axes x and y, and h (implicitly equal 
to 29) is the rotation angle about the third (z) axis. The 
superscript -I- or — on & tells how the angle varies as 
we move from one layer perpendicular to z to the next 
layer. Here, since there is only one layer, this super- 
script is meaningless. More generally, for RP systems 
we have two slabs to consider, there would be two sets 
of Glazer symbols, one for each slab. However, we will 
give a simpler symbol which applies when the star of the 
wave vector is specified. (The Glazer symbol implicitly 
specifies the wave vector by the array of superscripts.) 

For a rotation about the tetragonal z axis, the situ- 
ation is that shown in Fig. 01 There one sees that a 
vertex common to two adjacent octahedra would, if the 
vertices were considered to be separate vertices for the 
two octahedra, become two closely, but distinct, sepa- 
rated points. In order to recover the common vertex, the 
two points would have to coalesce, which would require 
octahedral distortions of Aa;/2 and Ay/2. However, it is 
possible for the lattice to relax, so that this mode would 
take place without any distortion of the octahedra. This 
relaxation involves microscopic strains along the x and y 
axes, to account for the A displacements. In the presence 
of microscopic strains Cxx and e^y, one has for RP214 

^x,n = a[29'^ + Cxx] , 
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FIG. 4: (Color online) The octahedral constraint for in- 
terlocking ^-rotations about the tetragonal z axis. Here 
Ax — Ay = 2a6^, where Ac is the mismatch that has to 
be absorbed either by a distortion of the octahedron or by a 
macroscopic strain when the octahedron is undistorted. 



a[29l 



(5) 



where Ac^„ is the value of A^ for the nth slab and 9n 
is the value of 9 for the nth of the 2 slabs in the RP214 
system. Therefore for RP214 the free energy per octahe- 
dron for the octahedrally constrained ^-rotated structure 
contains the term 
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where we have introduced the expansion parameter A 
which is the ratio of the stiffness coefficient for distorting 
an octahedron to other stiffness coefficients of the lattice 
and ce, and below c^, are constants of order unity. Most 
of our results will carried only to leading order in 1/A. 
Corrections higher order in 1/A are discussed in the Ap- 
pendix. Note that at order 1/A it is inevitable that the 
octahedra will in fact be distorted. 

We need to generalize Eq. ([6|) to allow for the rotation 
of octahedra about the tetragonal x and y axes. In the 
seminal work of Ref. [l3 we are reminded that the group 
of rotations is nonabelian, i. e. rotations about different 
axes do not commute with one another. Here we discuss 
a simple nonlinear treatment of small rotations. We de- 
fine the orientation of an octahedron by three variables, 
9, (j)x, and 4>y which correspond to small rotations about 
the tetragonal axes.[26j These three variables correspond 
to the three Euler angles needed to specify the orientation 
of a rigid body when the rotation from the undistorted 
state is small. This situation is somewhat similar to the 



spin wave expansion in which one introduces transverse 
spin components to describe the rotation of a three di- 
mensional spin. A simple way to deal with this situation 
is to express the orientation of the vectors from the cen- 
ter of the octahedron to the equatorial vertices in terms 
of their transverse displacements as 



ri = 
r4 = 



(^\/l/4-y?-z2,yi,Zi^ , 



-A/l/4-2/|-2|,y3,Z3 



X4, - 



1/4 -xl- zl,Z4 , 



(7) 



where ri, r2, rs, and r4 are the rotated positions of 
the vertices whose respective original locations were 
(1/2,0,0), (0,1/2,0), (-1/2,0,0), and (0,-1/2,0). 
Clearly, to retain the octahedral shape (with the center of 
mass fixed) we require that = — ri and — — r2. We 
wish to incorporate the octahedral constraint to leading 
order in the transverse displacements. This constraint 
leads to 



{Zx ± Z2)^8) 



so that 
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This gives 







4+yi\l^-X2-4 + Z1Z2 ■ (9) 
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To make contact with the body of the paper replace zi 
by (j)^ and Z2 by (j>y, to get 
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+ ZiZ2 + Oii^P) 



where q is one or more of the variables. Thus 

^^±^ = ^z,Z2 + X2{yl + zf)+y,{xl + zl) + Om) 

We transform from the variables X2 and yi to 9 and 69: 

-9 + 59; yi=9 + S9, (12) 
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so that 



59 = -ziZ2 + [-9 
+ [9 + 59][{-l 



^ 59][{9 + 59f + zl] 
+ 59f + z|] + O(g^) , (13) 



which gives 
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One sees that in terms of the variables 9^ clull yjy^ 
actual positions of the vertices are given by power series 
expansion in these variables, the lowest term of which 
identifies these variables directly with the corresponding 
displacements. (Of course, this expansion is only useful 
if the variables are small.) We can use these variables 
to analyze the symmetry of the free energy. In calcu- 
lating scattering cross sections one must, of course, use 
the actual positions of the ions given by their nonlinear 
expansion, such as that in Eq. or in Fig. [SJ below. 
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FIG. 5: (Color online) The displacements of the vertices for 
the star of X including nonlinear contributions. The displace- 
ments of the lower left octahedron at the origin are given in 
terms of the expansion parameters of Eq. (|15p . The param- 
eters at each vertex are shown. The 9 displacements are only 
given in the formulas. Because the wave vectors of the star 
of X are (1/2, 1/2,0) and (1/2, -1/2,0), the displacement of 
vertex C is related to that of vertex A by changing the signs 
of all the order parameters. The displacement of vertex F is 
obtained similarly from vertex D. The displacement of vertex 
B is obtained from that of vertex C by inversion about (1, 0, 0) 
and that of vertex E from that of vertex F by inversion about 
(0, 1,0). The origin is indicated by the filled magneta circle. 



In Fig. [S] we identify the displacements of the ver- 
tices for the star of X. From this figure we obtain the 
mismatch of the vertex when considered to be two inde- 
pendent vertices of adjacent octahedra to be 



xb - xa = 2(j)l 



2r 



dux 
dx 



ys - yA 



2(j)x 



dUy 

dx 
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XE - XD 



duj. 

2(px4>y + -TT- : 

dy 



(16) 



with corrections at order q^. Here u{v) is the displace- 
ment field whose derivatives give rise to the strain tensor. 
Thus, for the star of X we have 
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where e^y = [dux/dy + duy/dx]/2. Note that it is neces- 
sary to invoke a macroscopic strain to obtain a free en- 
ergy of distortion which does not involve distorting the 
octahedra. [mi Coupling the 9 or (f) variables to a trans- 
lational phonon will not bring all the octahedra along an 
axis closer together. 

One might argue that "Constraining the form of the 
free energy departs ... from the accepted way of using 
symmetry in this theory. Any author is certainly free to 
postulate a modified free energy and derive consequences 
but the general appeal of this approach is then limited 
...."[131 The reason this objection is invalid is that our 
treatment is predicated on the fact that these RP systems 
consist of rigid octahedra. (This assumption of rigidity 
has been accepted by the research community for sev- 
eral decades and is supported by recent first principles 
calculations. [6, 28- 30j ) One can imagine raising the tem- 
perature sufficiently or reducing the internal force con- 
stants so as to violate our assumption that the quartic 
potential due to intraoctahedral interactions dominates 
the quadratic terms in the Landau expansion. We refer 
to this limit as the limit of "octahedral melting." For the 
RP perovskites, this limit is clearly irrelevant in prac- 
ticality. But in perovskites, such as RP214, it is the 
geometry of the metal-oxygen bonds that leads to the 
rigidity of the octahedron. To ignore this physics and 
rely solely on symmetry (as implied by Ref. '27^ is not 
sensible. To summarize: if it is legitimate to consider the 
system as consisting of rigid octahedra (as is the case for 
the RP perovskites), then the elastic energy quartic in 
the ionic displacements is dominated by coupling terms 
which arise from the distortion of individual octahedra. 
Note that the octahedra do not need to be infinitely rigid 
for our argument to be valid. They only need to be rigid 
enough that the parameters of the quartic potential are 
not very different from those for rigid octahedra. 



III. RP214 STRUCTURES 
A. The star of X 

As mentioned, we will develop a Landau expansion 
for RP214 structures associated with the star of the 
wave vector, X, which includes Xi = (1/2,1/2,0) and 
X2 = (1/2, —1/2,0) providing that the octahedra rotate 
as constrained by their shared vertex. (See Fig. [5]). In 
the z = plane both Xi and X2 each imply that all 
angular variables alternate in sign as one moves between 
nearest neighbors. This fact fixes the values of all the 
0's and 6''s in the z — plane in terms of the values 
for the octahedron labeled A in Fig. [51 If we had only 
the wave vector Xi, then the variables for octahedron B 
would be the negatives of those for octahedron A and the 
variables of octahedron C would be identical to those for 
octahedron A. If we had only the wave vector X2 then the 
variables for octahedra B and C would be reversed from 
what they were for wave vector Xi. Thus, if we have 
a linear combination of the two wave vectors, the orien- 
tational state for the plane z — 1/2 is characterized by 
assigning arbitrary values to the variables of octahedron 
B relative to which the values of all the other variables in 
that plane are fixed. Thus Fig. IH] gives the most general 
structure associated with the star of X. 

Therefore we write the elastic free energy for the 
RP214 structure of Fig. [Blfor the star of X as 

2 

F{ct>x,kAy.k,ek) = cea^A^ [(202^^ + 2f?2 + e^^)2 

fc=i 

+ {2<Plk + 291 + fyy) 
2 

-l-c^a^ A ^ {2(l>x,k<f>y,k + ixy)'^ 

k=l 

+F2 +Fi+F,+ F2,e , (18) 

where the first terms come from Eq. p7)) and last line 
contains terms of order A°. Here F2 {F4) is the free energy 
quadratic (quartic) in the angles 6 and 0, -F^ is the strain 
free energy at quadratic order, and i^2,e is the rotation- 
strain coupling which is quadratic in the rotation vari- 
ables and linear in the strains. Here we ignore quadratic 
terms involving modes which distort the octahedra, since 
they will not be activated. Also, as we shall see, terms 
coupling (j) and 9 variables such as ^9l{(j)1,, + (j)y i,) do not 
affect the results because 9 and cj) are not simultaneously 
critical. 

The free energy has to be invariant under all the sym- 
metry operations of the "vacuum," which, in this case, 
is the parent tetragonal structure. Accordingly, in Table 
U] we give the effect of symmetry operations on the vari- 
ables appearing in Eq. (IT^ . Here and below, because 
of the octahedral constraint quartic terms of the form 
9i9h 4>'l,i4>l,2 + 4>i,24'li, and (l>x,i(l>y,i<l>x,2<l>y,2 which are 
allowed by symmetry (see Table I) do not appear at order 
A. (But, of course, they are present at order A°. Only 
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FIG. 6: (Color online) The structure of corner-sharing octahe- 
dra for the star of X = (1/2, 1/2,0) and P = (1/2, 1/2, 1/2). 
The solid squares are the cross sections of octahedra in the 
z = plane and the dashed squares are in the z = 1/2 plane. 
For clarity the octahedra are slightly separated instead of 
sharing vertices. Here (px means that the +x vertex moves 
up by an amount (j)x and the —x vertex moves down by an 
amount (px and similarly for (py- Also 26 is the angle of rota- 
tion about the z axis. Here and below Q denotes —Q. For X 
the structure is invariant under z — > z + 1. For P the variables 
change sign under z ^ z + 1. The octahedral labels A, B, and 
C are needed for the discussion in the text. 



quartic terms which arise from the octahedral constraint 
at a single vertex can appear at order A.) Using Table U 
we see that the quadratic terms which are invariant un- 
der the symmetry operations which leave the reference 
tetragonal structure invariant are 
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in the Voigt notation[3l| where 1 = {x,x) 
etc. Similarly, we use Table U to write 
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We will deal with F^ when it is needed. We will give 
analysis of the above free energy which neglects fluctua- 
tions. Our model has some resemblance to that of Bean 
and Rodbell 32] except that the form of the free energy 
does not drive the system to a first order phase transi- 
tion, at least within mean field theory. An interesting 
problem would be to give a renormalization group anal- 
ysis like that of Bergmann and Halperinjsl] to elucidate 
the effects of orientation-strain coupling on the structural 
transitions. 



TABLE I: Effect of symmetry operations on the variables of 
the stars of X = (1/2, 1/2, 0) and P = (1/2, 1/2, 1/2) (shown 
in Fig. [6]) and on the strain variables. Here TZi is a four- 
fold rotation about the tetragonal z-ajds passing through the 
origin, nid and are mirrors that take x into y and z into 
—z, respectively and T is the translation (1/2,1/2,1/2). These 
variables are odd under the translations = (1,0,0) and 
Ty = (0, 1, 0). For the star of X, ^ = 1 and for the star of P, 
^ = — 1. Note that spatial inversion I is implicitly included 
because T — 7?.^m^.[34l| In the last line = e^x — f-yy 
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The structural phase transitions which we are investi- 
gating arise when one of the channels becomes unstable, 
i. e. when 7 or a — passes through zero. (As in Ref. 
[TtI . we reject multicritical points where more than one 
channel simultaneously becomes unstable.) 



1. 9 distortion 
^p) For instance, when only 7 becomes negative, then 
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where F4 {e) are quartic terms in e of order A" . To leading 
order in 1/A, when F{e) is minimized, one finds that 
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FIG. 7: (Color online) As Fig. [G] The structure of corner- 
sharing octahedra obtained by 9 rotations either for X = 
(1/2,1/2,0) (top) and P = (1/2,1/2,1/2) (bottom). The 
arrows indicate the displacement of the oxygens in the equa- 
torial plane. Here and in the figures below the original tetrag- 
onal axes are labeled by lower case letters and those of the 
distorted structure are labeled by capital letters. The ma- 
genta dot represents the tetragonal origin. In the top (bot- 
tom) panel X = (0, 0, l)t {Z = (0, 0, 2)t), where the subscript 
indicates that components are taken in the original tetragonal 
system. In the top (bottom) panel the distortion is unchanged 
(changes sign) for z ^ z -\- 1. In the top (bottom) panel the 
new origin is at 2: = (z = 3/4). 



(Corrections to this result at leading order in 1/A are 
studied in Appendix A. In particular, there are no such 
corrections to the relation d1 = 62 ) As we discuss in a 
moment, all choices of signs for the 6''s in Eq. (|24|) lead 
to equivalent structures, one of which is shown in Fig. 

Ell. 

Now we identify the space group for the above 9 
distortion. The X- structure has generators (s^l {X ± 
1/2,Y + II2,Z), (X,y,Z+l), (X,F,Z), (X,F,Z), and 
{X + 1/2,Y,Z + 1/2). In determining the space group 
from these generators, j36| it is useful to realize that the 
structures we find must form a subset of those listed in 



Ref. [13. We thus identify the space group of the struc- 
ture of Fig. [7^ as or Cmca (#64). Cmca (64), one of 
the three ^-dependent structures of irrep for the star 
of X which are listed in Ref. However, we do not al- 
low the other two structures of Ref. [l7|, the first of which 
(Dt^ or Pbam=#55) has, according to Table III of Ref. 
[itI . |6'i| 7^ 1 6^2 1, with 6*1 6*2 7^ and the second of which 
(D|^ or P4/mbm=#127) has, according to Table III of 
Ref. [13, one sublattice distorted and the other not dis- 
torted, so that 01^2 = (see Fig. [3]). The problem with 
these structures is that to avoid distorting the octahedra 
we had to invoke a uniform strain which only relieves the 
distortions in the two slabs when the distortions in the 
two slabs are the same. Thus when 6\ ^ O"^, there is an 
unavoidable distortion energy of the octahedra of order 
A. 

Note also that the term in proportional to txyOi02 
in combination with F^^ leads to 

f-xy -ai0i92/cii , (25) 

This distortion is consistent with orthorhombic symme- 
try and with the orientation of the orthorhombic axes 
shown in Fig. [T] Rotating the crystal by 90° {Ti-i) 
changes the sign of 9i92 and thus changes the sign of 
txy, as one would expect. 

In comparison to other ordering transitions we can 
make an analogy between the order parameters which 
govern the distortion from the parent tetragonal phase 
and the order parameters in, say, a magnetic system. 
In this formulation the distortion of the parent lattice 
in pcrovskites is analogous to the development of long- 
range magnetic order. Having a distortion only within 
one sublattice of the RP system is thus analogous to 
having magnetic order only on some sublattices. Al- 
though one can have ordered systems which have some 
disordered components, they differ from the present case. 
For instance, the orientational phase II of solid methane 
(CD4) consists of a unit cell having six orientationally or- 
dered methane molecules and two completely disordered 
molecules. [37} In that case, the site symmetry of the dis- 
ordered molecules is high enough that the effective field 
from the ordered molecules vanishes. Furthermore, the 
interaction energy E between two disordered molecules 
is much less than kT, even at the lowest temperature 
at which phase II exists, so that they do not coopera- 
tively order. Superficially, the situation here is similar to 
that for solid methane. For instance, suppose one builds 
up the distorted structure plane by plane for wave vec- 
tors in the star of X or P. The first plane would have 
9 = 9o, say. Moving to the second plane we note a frus- 
tration due to the four-fold rotation, TZ^ which implies 
that the energy is invariant against changing the sign of 
92- (Table H] indicates that TZ4 leaves 9i fixed but takes 
6*2 into —02.) If one considers the third plane, there is 
no analogous frustration. One will have 9^ = ±6*1, the 
sign depending on the weak interaction between octahe- 
dra in plane #1 and those in plane #3. More generally, 
0n+2 — o'9n, where a can be either -1-1 or —1. As in quasi- 
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two dimensional systems, as the temperature at which 
the distortion becomes unstable is approached, two di- 
mensional correlations will become significant and then 
even a weak coupling in the third dimension will lead to 
three dimensional long-range order at a single phase tran- 
sition. This situation is reminiscent of the decoup ling of 
magnetic sublattices in the bcc antiferromagnet. |38ll39l | 
The point is that if the distortion order parameter be- 
comes nonzero in, say, the even numbered planes, the 
mechanism that led to this order would also apply to 
the odd numbered planes, which would then also distort 
at the same time. The distorting of both sublattices of 
octahedra could only be avoided if simultaneous distor- 
tions were strongly disfavored by the form of the quartic 
potential (i. e. if w of Eq. Q were positive.) This pos- 
sibility seems unlikely and indeed our analytic treatment 
of the octahedral constraint indicates that this scenario 
does not occur for large A. 



a) bi 



18 

D2h= Cmca = 64 




2. (f) distortions 

Now drop the variables, so thatji^ the (^-dependent 
free energy for the star of X is 

2 

k=l 

[{a - /?)/2] - 0,,2)2 + ((/.,,2 - 

+ [{a + l3)/2] + ^y^2f + (0..,2 + 

2 

+C4,0?\ ^ i2(f>x,k4'y,k + £xyf + ^4(0) , (26) 
fe=l 

where Fi((j)) are the terms quartic in (f> which are pro- 
portional to A''. When F{(j)) is minimized for large A, we 
get 

<^xx = -S^^^i = -'^4>'i,2 ^ 4'x,2 = ^4>x,l , (27) 

eyy = -20y^i = -20^ 2 =^ <t>y.2 = ±'/>y,i , (28) 
and 

e^y = -20j;,i(/)y_i = -24>x^24'y,2 ■ (29) 

There are four possible directions of the ordering vector 
= [(l)x,ii (t>y,iT<t'x,2, (t>y,2] depending on whether or not 
a — f3 becomes critical (zero) before a + /3 and the choice 
of signs in Eqs. (I?7l) and (1^ . When a — (3 is critical (so 
that 4'x.i + 4'y,2 = 4'x,2 + 4'y.i = 0), then is proportional 
to either ai = [UTT] or 61 = [UlT]. li a + (3 is critical (so 
that 4'x,i—<Py,2 — <Px,2 — 4'y,i =0)1 then ^1/ is proportional 
to either a2 = [1111] or &2 = [HH]- In each case, the two 
choices are equivalent: [4l| 7?.46„ = a„. Figure 0jhows a 
representative of these solutions for each case. 42] 

Finally, we identify the space groups of the structures 
of Fig. E The generators of 61 are {X ± 1/2,Y + 



b) a2 ^1 D2h= Cccm = 66 




FIG. 8: (Color online) As Fig. |6] for the star of X (with 
invariance under z z + 1). x, y, and z are the tetragonal 
axes and and X,Y, and Z are the conventional lattice vectors 
after distortion. The filled magenta circle is the tetragonal 
origin. * = [(px,i,(l>y,i,(l>a:,2,<j}y,2] = [iTlT] for a) and [1111] 
for b).[42] The irreps are given in Table Hill The new origin 
has z = for a) and z — 1/4 for b). The third new axis vector 
is [001] t in terms of the original tetragonal coordinates. 

l/2,Z),(X,r,Z + 1), (X,F,Z), {X,Y,Z), and (X + 
l/2,r, 1/2 + Z) _and_those_of 02 are {X±l/2,Y+l/2,Z), 
{X,Y,Z+1), {X,Y,Z), {X,Y,Z), and {X,Y,l/2 + Z). 
From these generators we identify the space groups as in- 
dicated in Fig. m From Eq. (|2T|) we see from the term in 
as that exx = ^yy for both these structures. The term in 
og indicates that exy ^ 0, consistent with the orientation 
of the coordinate axes shown in Fig. El 

Note that, in comparison to Ref. llZl, our formula- 
tion does not allow the structures of space groups Pccn 
(Di° = #56) and Pnnn (D^^ = #48). From Table III of 
Ref. [l^one sees that Pccn has 4>x.\ = 4'v,2 = a and ^y^\ = 



8 



02;. 2 = b, with a 7^ 6 and Pnnn has (j)x,i = ~4'y,2 = a and 
<l)y^i = —(j)x,2 = with a ^ b. As before, to reheve the 
distortion of the octahedra via strains imphes that the 
order parameters have the same magnitude on both sub- 
lattices. Similarly we do not allow space groups P42/ncm 
(D^s = # 138) and P42/nnm (D^^ = # 134) which have 
(in one setting) (t)xA = 0y,2 = and \4ix,2\ = \4'y,i\- Our 
analysis of the space groups listed in Ref. [l3 is summa- 
rized in Tables HI] and inil 

TABLE II: Space groups of Ref. [Ij for the stars if X, N, 
and P which we do not allow for RP214's. For the column 
headings "S " stands for the Schoenflies symbol, "H-M" stands 
for the short Hermann-Maugin symbol as given in Ref. 8, and 
# is the number of the space group in Ref. g. For otherwise 
identical space groups, the footnotes give the basis vectors of 
the unit cell in terms of the original tetragonal coordinates. 
The irrep labels are from Ref. 1 171 . 



Irrep 


S 


# 


H-M 


S 


# 


H-M 






127 


P4/mbm 


'-'2h 


55 


Pbam 






56 


Pccn 


p)16 


138 


P42/ncm 






48 


Pnnn 


T^12 


134 


P42/nnm 
PT*** 


N+ 




10 


P2/m 


c\ 


2 


N+ 


T^2S 
^2h 


74 


Imma* 


T-,25 
^2h 


71 


Immm 


N+ 


'^2h 


15 


C2/c 


p3 
"-'2h 


12 


C2/m1" 


N+ 


p3 
'^2h 


12 


C2/mtt 


p)19 

^Ah 


141 


I4i/amd 


N+ 


'->2h 


65 


Cmmm 


'-'Ah 


139 


I4/mmm 


P4 




140 


I4/mmm 


^2d 


122 


I42d 


P5 




74 


Imma** 


^2h 


69 


Fmmm 


P5 


Si 


82 


14 


(^22 
'-^2i> 


46 


Ima2 


P5 


'^2v 


43 


Fdd2 


^2v 


42 


Fmm2 


P5 




24 


I2i2i2i 


Di 


23 


1222 


P5 




22 


F222 


^2h 


15 


C2/c 


P5 


'^2h 


12 


C2/m"t 




5 


C2 



*: (200), (020), (002). **: (002), (110), (110). 
_ _ (Tll)_, (iii, (llT). f: (002), (220), (Til), 
tf: (202), (020), (0,0,2). \\\: (020), (002), (110). 



B. The star of N 

Similarly, we construct the most general structure for 
the star ofN = (1/2, 0, 1/2) whose wave vectors are Ni = 
(1/2,0,1/2), N2 = (-1/2,0,1/2), N3 = (0,1/2,1/2), 
and N4 = (0, -1/2, 1/2). To do this we use Fig. El Con- 
sider first the situation in the 2 = plane. Note that 4>x 
has to alternate in sign as we move along x. This means 
that is associated with a linear combination of Ni 
and N2 distortions. Since N\^y = iV2,j, = 0, we see that 
(\)x must be independent of y. Similar reasoning indicates 
that 0y alternates along y and is therefore associated with 
a linear combination of N3 and N4 distortions. There- 
fore is independent of x. These wave vectors do not 
support nonzero values of d. We have thereby fixed all 
the values of the variables in the = plane in terms of 



TABLE HI: As Table El Space groups of Ref. |T3 for the 
stars of X, N, and P which we do allow for RP214's. Under 
"Var's" we give the variables active in the mode and under 
"Fig" we give the number of the illustrative figure. 



Irrep 


Var's 


s 


# H-M 


Fig 


x^ 


e 


^2h 


64 Cmca 


7a 


xj 




1^18 
^2h 


64 Cmca 


8a 




(t>x,(f>y 


T-)20 


66 Cccm 


8b 


N+ 


4>y 


p3 


12 C2/mt 


10b 


N+ 


<i)x,(py 


p3 


12 C2/m" 


10c 


N+ 


<t>x, <t>y 


c} 


2 PT" 


10a 


P4 


9 




142 I4i/acd 


7b 


P5 


<t>x,(py 


T^2(i 

^2h 


72 Ibam 


Ua 


P5 


<i)x,<t)y 


f~*6 
^2h 


15 C2/c 


11c 


P5 


4>y 


^2h 


70 Fddd 


lib 



f: (Oil), (100), (011)_. 
tt: (002), (220), (i|i). 
* The wave vector may be close to, but not exactly at, the 
star of P. 
**: (Til), (iTl), (llT). 




FIG. 9: (Color online) As Fig. [6l The structure of corner- 
sharing octahedra for the star of N. The variables change 
sign under z ^ z + 1. 



those of octahedron A. Now consider the situation in the 
z = 1/2 plane. Suppose we have a linear combination 
of Ni and N2 which gives rise to (j)x,i for octahedron A. 
If we had only Ni, then the variables for octahedron B 
would be —4>x,i and —(t>y,i, whereas if we had N2, then 
these variables would be (px.i and 4'y.i- As for the case 
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TABLE IV: As TableU] but for the RP327 variables of the star 
of N = (1/2,0, 1/2) (shown in Fig. HJ. T = (1/2, 1/2, 1/2), 
T:, = (1,0,0), and Ty = (0, 1,0). All variables are odd under 
= (0,0,1). 





7^4 






r 




Ty 


4>x,i 




<Py,i 


— <t>x,l 


<t>x,2 






't>V,l 






-<l>y,i 


4'y,2 


't'y,l 


-4'y,l 


4>x,2 


4'y,2 


't>y,2 


4>x,2 


4>x,i 


— 4>x,2 


4>x,2 


<t>y,2 


4>x,2 


<t>x,2 


0y,2 




0y,2 





of the star of X, we conclude that all the variables in the 
second layer are fixed in terms of the arbitrary variables 
of octahedron B, so that Fig. ^ characterizes the most 
general structure arising from the star of N. The effect of 
symmetry operations on these variables is given in Table 

llYl 

We need to establish the analog of Fig. [5] for the star 
of N. First of all, here we can set 9 = 0. Also, note that 
for sites B and C, the sign of (f)y is reversed for the star 
of N in comparison to that for the star of X. For site A, 
r^i for the star of N is as in Fig. [5] However, now 



tb = {1/2 + (j)l,-(l)^(j)y,(j}x) 



(30) 



Therefore ub ^ VA ^ zb — za = and xb — xa = 20^. 
Likewise, for site D Fig. [5] applies equally for the star of 
N. However, for sites E and F the sign of (j)x is reversed 
for the star of N from what it was for the star of X. So 



Ye 



(31) 



20? 



so that xe — xd = ze ^ zb = and i/e — yo - --ry 

We now write the Landau expansion for the star of N. 
Taking account of the symmetries of Table IIVI and the 
preceding discussion, we obtain the relevant terms in the 
free energy for the star of N to be 



F, 



N 



cga A 
+ H 



2 

x.l 



2 

x,2 



-f[<r 



^VVJ 



2e 



where 



u 



- 4 
+v 



-X \ 



i'll + > 



+ 



'^,20y,2) ^ 



^2) 



^lA<fi,2 



Fiim 



(33) 



To leading order in 1/A the minima of F^ occur for 



Ma 



y-2 



'Jy.2 



•x.2 — ± 
± 



>xA 



(34) 



The quartic terms distinguish between these solutions. 
If (/)^ = </)^ ]^ = 0^ 2 similarly for <^y , then we have 



F4 



2{v-w + x)c\>ic\>'y + [ii; + {ul2mi + 0^)^ (35) 



This indicates that we can have a structural phase tran- 
sition into three classes of states. We can have a contin- 
uous phase transition into states of class A with = 0^ 
if (ti — w -I- x) is negative or into a state of class B with 
'\>x'i>y = \i (y — w -\- x) is positive. In addition, wc can 
have a phase transition to a state of class C in which 
(\)x and (^y do not assume special values if the higher or- 
der terms cause the transition to be discontinuous. The 
next step is to determine which of these solutions are in- 
equivalent jL e. which are not related by a symmetry 
operation. 41 

We first show that all solutions of class A are equivalent 
to one another. Using the results given in Table IIVI we 
have that 

(l + 7^4 + 7^^ + 7^|) [1111] = [{mj^Ih] , (36) 

where {/ii^} indicates the set of v values, i. e. {^i^} = 
11 + IT + Tl +TT. Then 

T(l + 7^4-f 7^^ + 7^;^) [nil] = [n{/xi/}] = $ ,(37) 



and finally 



^1 + 7^4 + + Te^) $ = [{pr}W] 



(38) 



where (apart from an arbitrary amplitude) the right- 
hand side of this equation includes all vectors of class A. 
In class A, i/)^ — ^"i/i so that exx — ^yy and we take [TlTT] 
as its representative. The term in Eq. (|2ip proportional 
to ag indicates that Cxy ^ for this structure. 

Next we consider solutions of class B. From Table ITVl 
note that 



and 



^4][1 



7^4<I> 



r,][ioio] 



E ([/iOi^O] -I- [OAiOi^] 



$ . (39) 



(40) 



The right-hand side includes all vectors of class B. We 
may take [OlOT] as the representative of class B. 

Finally we consider solutions of class C, which are of 
the form [a;, y, ±x, ±?/], where ^ \y\ and both are 
nonzero. Using Table HVl we write 



'mdTlA\x,y,x,y\ = {xyxy\ 



(41) 



so that 

[1 + m. 



\^ + mdTl^{xyxy\= ^ ^ (crix,CT2y,a;,?/](.42) 

CTi— ±1 cr2— ±1 



From this we conclude that all vectors of class C are 
equivalent to one another, and we take their representa- 
tive to be [S'yxy]. 

Thus, in all, we have the three allowed space groups 
from the star of N shown in Figs. [TU] and [Tl] [xyaJy], 
[OlOT], and [TlTT]. As before, to determine the space 
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groups of the structures we use Figs. [TOl and [TT] to write 
down their generators. For [xyxy] the generators are 
{X + 1,Y,Z), {X,Y + l,Z), {X,Y,Z + l), and (X,F,Z), 
which is G\ (PI = # 2). For both [OlOT] and [TlTT], 
the generators are (X,Y^_+ 1), (X + l/2,r_+ 1/2, Z), 
(X + l/2,r-l/2,Z), and which is 

(C2/m = #12). Ahhough these two structures be- 
long to the same space group, they are different because 
their unit cehs are different (see Figs. [TOl and fTTj) . 

In Table IIIII (ITU we hst the structures which are (not) 
aUowed. We only allow one of the eight structures for the 
irrep listed in Ref. [TtI which occur via a discontinuous 
transition. In addition, Ref. [O lists two space groups 
Cmmm (65) for which * = [1000] and I4/mmm(139) 
for which ^ = [1100]. Both these are inconsistent with 
the fourth order terms arising from the rigid octahedral 
constraint. [I. e. they did not arise from Eq. ([5^ ]. 
They are also counterintuitive in that they both describe 
states in which nonzero order parameters appear only on 
alternate planes. [Look back at the discussion below Eq. 

mi 

C. The star of P 

The star of P consists of the vectors Pi = 
(1/2,1/2,1/2) and P2 = (1/2,-1/2,1/2). The possi- 
ble structures are identical to those for the star of X, 
except that the variables change sign under z ^ z -\- 1 
as indicated in Fig. |6l The transformation properties of 
these variables are given in Table HI The quartic terms 
are the same as for the star of X. However, the quadratic 
terms differ because of the ^ factors that appear in Table 
m The octahedral constraint assumes the same form as 
for the star of X because the stars of X and P only differ 
in how the layers are stacked. Therefore for the star of 
P we have 

2 

F{4>x,k,4>v,k,ek) = cea^A^ [(202 ^ + 20^ + e,,)^ 

fc=i 

2 

-|-C0a^A^ {2(px,k(f>y,k + <^xvf 
fc=i 

+^^2 +Fi+F,+ F2,e . (43) 

1. 6 distortions 

We analyze these as before. In the channel where 7 
passes through zero, we have 

2 

F{e) = cea'xY, [(2^?^ + ^..f + {K + ^yyf 

k=l 

'Im&i + di] + FM + + , (44) 



c\ = 2 \y pi 




FIG. 10: (Color online) As Fig. [8]for the star of N (with sign 
change under z — !> z+1) for a) (class C) \xyxy\ and b) (class B) 
[0101]. X, y and z are the tetragonal axes, and X, Y , and Z are 
[(Til), (1/2, -1/2, 1/2), (UT)] in a) and [(OTl), (100), (Oil)] in 
b). All the new origins are at 2 = 0. 



SO that \9i\ = \d2\- This structure has the same degen- 
eracy associated with the relative phase of even and odd 
layers that we saw for the previous structures (at the 
star of X). Because — 1/2, the only 0-dependent 
structure has 0n+2 = —On- The other two structures 
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c) \y 




FIG. 11: (Color online) Continuation of Fig. 1101 for c) (class 
A) [TlTT]. X, Y, and Z are [(002),(220)(-l/2, 1/2, -1/2)]. 
The new origins is at 2: = 0. 



listed in Ref. [TtI which have \9i\ ^ \92 \ are not admissible 
due to the form of the term proportional to A in Eq. (H^ . 
The allowed P-structure, shown in Fig. [Tja has generators 
{X - 1/2, r + 1/2, Z+ 1/2), (X_+_l/2, 1/2, Z+ 1/2), 
{X + 1/2, Y + 1/2, Z - 1/2), {X, r, Z), {X, Y,Z + 1/2), 
and (y + 1/4, X + 3/4, Z + 1/4), which is space group 
or I4i/acd (# 142). 



1;. (/) distortions 

Now consider the (/)-dependent solutions which are as- 
sociated with irrep P5 according to Ref. 1 7. Note that 
all the subgroups from this irrep listed there and in Ref. 
[T9I do not satisfy the Lifshitz condition. What this means 
is that the quadratic instability occurs at a wave vector 
that is not fixed by symmetry [43[ (so that it can not be 
assumed to be at the star of P). Accordingly, the wave 
vector can only be at the star of P if the transition is 
discontinuous. 

However, we can determine which structures with wave 
vectors either at or near the star of P can be condensed. 
To do this, we simply ignore the Lifshitz criterion in our 
analysis of the free energy. For this purpose we consider 
the (/)-dependent free energy which is 

2 

F = cea^Xj2 [(20^ + 6,,)' 

k=l 

2 

k=l 



+F4 + F, + F2,, . (45) 

and we consider the phase transition which occurs when 
a passes through zero. For large A we must have 

e^^ = -202 1 = -2(/)2 2 

eyy = -20^,1 = -20y,2 

txy = -'^4'x,l4>y,l = -201,, 2(^^,2 ■ (46) 

Thus we have three classes of solutions: in class A we 
have [1,0, ±1,0] and [0, 1,0, ±1]. In class B we have 
[llcrtr] and [IIctct] with a — I, class C is like B, except 
that < (7 < 1. These three classes are associated with 
the different minima of ^4 which assumes the same form 
as in Eq. ([33)) . The analysis to show that all members 
of same class are actually equivalent follows the previous 
arguments, so we will not repeat it here. The three solu- 
tions are shown in Fig. [121 We take the representatives 
to be [OlOT] for class A, [Till] for class B, and [xyxy] for 
class C. 

As before, we identify the space groups of these so- 
lutions by determining the generators of the represen- 
tatives shown in Fig. [T^l For a) ([Till]) the gen- 
erators arei44l iX,Y,Z), (X,F,Z + 1/2), (X,Y,Z), 
{X - 1/2, r-l- 1/2, Z-f 1/2), {X + 1/2, r- 1/2, Z+ 1/2), 
(X + 1/2,Y + 1/2, Z ~ 1/2) and thus the space group 
is Ibam {Dll =#72). For b) ([OlOT]) the genera- 
tors _a.T^(X + 1/4,F+ 1/4, Z), {X,Y + 1/4, Z + 1/4), 
{X, Y, Z), (X, Y + 1/2, Z + 1/2), (X + 1/2, Y, Z + 1/2), 
{X + 1/2, Y + 1/2, Z), and thus the space group is 
Fddd (Djl =#70). For c) {[xyxy]) the generators 
are (X,Y,Z + i), (X,Y,Z), {X + 1/2,Y + 1/2, Z), 
{X - 1/2, Y + 1/2, Z), {X,Y,Z +1) and thus the space 
group is C2/c (Cj/j =#15). When the Lifshitz insta- 
bility is resolved, these space groups give rise either to 
commensurate structures at the star of P via a first or- 
der transitions or to structures having incommensurate 
wave vectors near the star of P. 



IV. RP327 STRUCTURES 

In this section we perform the same analysis for the 
n = 2 RP systems A3B2C7 which we cah the RP327 
systems. As shown in Fig. [T] these systems consist of 
two slabs. Each slab consists of two layers which we 
label a and b (or 1 and 2). If we fix the (j) variables in 
layer a, the octahedral constraint fixes each variable 
in layer b to be the negative of its nearest neighbor in 
layer a. As a result each structure of the RP327 system 
is characterized by the same number of variables as its 
analog for the RP214 system. In contrast, since there is 
no such relation between the 9 variables of layers a and 
b, we introduce variables 9n^x to describe the rotation 
within the xth layer {x = a, b) of the nth slab (n = 1,2), 
as shown in Fig. 1141 The transformation properties of 
the variables are summarized in Table [Vl 
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a) D^^=Ibam = 72 [1111] 
2h 




FIG. 12: (Color online) As Fig. [8] for the star of P (variables 
change sign under z ^ z + 1). x, y, and z are the tetragonal 
axes and and X, Y, and Z are axes of the distorted structure. 
The new origins are in the z = 3/4 plane. The new out- 
of-plane lattice vector has magnitude 2 and is along z. The 
actual structures may involve an incommensurate wave vector 
near the star of P. 



A. The star of X 

1. 9 structures 

Using the symmetry operations of Table IVlv^fC find that 
the free energy of the 9 structures for the star of X as- 
sumes the form 



2 6 

F{e) = cga^Aj]^ [(2(?L + fxx)' + (2C + ew)'' 

fc=l a=a 



c) C,. = C2/c = 15 [xyxy] 




FIG. 13; (Color online) Continuation of Fig. [T2 The new 
origin is in the z = 1/4 plane. The new out-of- plane lattice 
vector has magnitude 2 and is along z. 



TABLE V: As Table H] for the stars of X and P, except that 
this table is for the variables of RP327 systems. Note that 
z = (about which rriz is taken) is midway between layers a 
and b. T is the translation (1/2, 1/2, 1/2). 





7^4 


rrid 


rriz 


T 


6l,a 


9l,a 


— 6l,a 


Oi,b 


62, a 


Ol,b 


Ol,b 


-6*1,6 


dl,a 


02,b 


02,a 


— d2,a 


— S2,a 






S2,b 


— d2,b 


— Q2,b 


^02,a 


^di.b 


4>x,l 


<f>V,l 




<t>x,l 


<t>x,2 


<Py,l 


— 4>x,l 


4>x,i 


't>y,l 


<t>y,2 


4>x,2 


-4>y,2 


<t>y,2 


£,4'x,2 


^4>x,i 


't>y,2 


4>x,2 


4>x,2 


i4'y,2 


^<l>y,l 



+ b [9ia0lb + 02a02b] + + F2, . 

For large A, the minima of this free energy as (a 
passes through zero occur for 



(47) 



-2ei 



-2^ 



26*26 



9lh , 92a — —-rr\^2b 



l^pro , ^.a- i^i^zo. (48) 

Thus, as shown in Fig. 1151 we have two possible dis- 
torted structures, depending on the sign of b: 

Ol^a = -Ol.b = e02,a = -C^2,6 , & > (49) 



l,a — 6*1, b — C^^2,a — ^0 



2,b 



b<0 



(50) 



fc=i 



where ^ = ±1. The fact that ^ can have either sign 
indicates the decoupling of even and odd numbered sub- 
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z = 0.6 



z= 0.4 



/(£) ® 


® / 











z = -0.1 



s 
1 
a 
b 



FIG. 14: (Color online) Values of the rotation (9) variables 
in different 2-planes for the most general such distortion of 
RP327's for the stars of X and P. The (red) dashed line is 
the axis about which the fourfold rotation TZ4 is taken and 
the (blue) dots show the points a; = y = in the z = 0.4 and 
z = 0.6 planes. Under TL^ an octahedron in slab #1 is taken 
either into itself or into an equivalent octahedron. Under TL^ 
an octahedron in slab #2 is taken into another whose rotation 
angle is of opposite sign. For the star of X, S„+2,a! = Qn,^ and 
for the star of P, Qri+2,x ~ '^Qn,x- 





FIG. 15: (Color online) The two 9 modes for a bilayer. Left: 
the two layers (a and h) rotate in phase as in Eq. (|50p . Right: 
the two layers rotate out of phase as in Eq. (|49|l . The out of 
phase rotation increases the energy by twisting the oxygen or- 
bitals but this is compensated by reducing the Coulomb inter- 
actions between octahedra. First principles calculations [30l| 
indicate that these modes differ only slightly in energy. 



cause undue difficulty in visualizing the effect of rotations 
about an axis perpendicular to the plane of the paper. 
However, operations such as reflection through the plane 
of the paper or inversion relative to the center of the octa- 
hedron introduce visual complications because these op- 
eration interchange inner and outer symbols. Experience 
indicates that for these operations one should use the 
results of Fig. [ITl supplemented, if need be, by a trans- 
lation. Operations such as O = (X, Y , Z) which involve 
taking Z into —Z are best expressed as O = X{X ,Y, Z) 
or O = vTLziX^Y , Z). The mirrors perpendicular to the 
page do not cause any confusion because for them outer 
symbols are taken into outer symbols, thus avoiding vi- 
sual complications. We now return to the discussion of 
the mode of Eq. (HHl) in Fig. [TBk . Using, if need be, the 
results of Fig. [ITl one sees that this structure has gen- 
erators (X_±J/2,y-f- 1/2,Z}, {X,Y,_Z + 1), (X,Y,Z), 
{X + l/2,Y,Z + l/2), and {X + 1/2,Y, Z), which there- 
#68 coming from irrep Xj^. 



fore is Ccca (Dj^J 



lattices which we mentioned in connection with the anal- 
ogous RP214 9 structures. 

We now determine the space groups corresponding to 
these two modes. In the mode of Eq. ([50)) the two lay- 
ers can be coalesced continuously into a single layer. So 
this structure has the same symmetry as the Cmca 9- 
structure resulting from the star of X and is obtained by 
replacing each layer by a bilayer with in phase rotations. 
See Fig. [7]a. The mode of Eq. ^ is shown in Fig. [TBb . 
We have to discuss the way we depict bilayer systems in 
our figures. Within each square (which represents the 
equatorial oxygen of an octahedron) the symbols closer 
to the corners of the square apply to the upper layer of 
the bilayer and the symbols closer to the center apply to 
the lower layer of the bilayer. This convention does not 



2. <j) structures 

Now we analyze the free energy for cj) distortions at 
the star of X. Here the effect of to^ has the opposite sign 
from the RP214 case. But since all terms are of even 
order, the symmetry of the free energy is the same as for 
RP214. So 

2 

F = cea'A^ [(2^2^-f + (2(/.2^ + eyy)' 

2 
fe=l 

+ [{a - /3)/2] [(0,,i - c^ya? + {^.a ~ 
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b) 




FIG. 16: (Color online) As Fig. [T] The 61-modes with 
On,a ~ —9n,b- The full (dashed) squares are equatorial sec- 
tions ol the bilayer centered at 2 = (z = 1/2). The outer 
full (inner dashed) arrows are the displacements of the equa- 
torial oxygens in the equatorial plane in the upper (lower) 
layer of the bilayer. In the upper (lower) panel the arrows are 
unchanged (reversed) under z z + 1. The new origin is at 
z = in the upper panel and at z = 1/4 in the lower panel. At 
right we show the z-coordinates of the full and dashed layers. 



where 4>a,n is the value of (pa for the top layer of the nth 
(n = 1,2) slab. 

The analysis parallels that for RP214 systems. Here we 
characterize the structures by the values of (f>n,a in the top 
layer of each bilayer. From the first line we conclude that 
'^Ki = '?^s2 ^iid (pyi = 0^2- There are two cases: the first 
is if a— /3 is critical and the second is if a -I- /3 is critical. In 
the first case (f)xi = ~4'v2 and 4>x2 = —<Pyi- In the second 
case 03;2 = 4>yi and = 4>y2- Then for the first case the 
possible ordering vectors are proportional to $ = [llTT] 
or [lllT]. These are equivalent structures and we take 
the second one as the representative for this case. For the 
second case the possible ordering vectors are proportional 
to $ = [1111] and [1111]. These are equivalent structures 
and we take [llllj as the representative for this case. Fig. 
[TH] shows these representatives. 




FIG. 17: (Color online) The effect of operations which take Z 
into ~Z. Here the distortion of a bilayer is represented by two 
sets of symbols within the square representing the equatorial 
oxygens. The outer set of symbols applies to the upper layer 
of the bilayer and the inner set to the lower layer of the bilayer. 
Inversion (X) is taken with respect to the (blue) filled circle as 
the origin, which for bilayers is in the plane midway between 
the upper and lower layer, is a mirror operation with 
respect to the plane perpendicular to the tetragonal z-axis 
which passes through the origin. 



Finally, we identify the space groups of the structures 
of Fig. [iHl Note that a -I- /3 critical corresponds to irrep 
and a — P critical corresponds to irrep in the no- 
tation Ref. ilOj. The generators of [iTlT] are (X±l/2, Y+ 
l/2,Z_UX,Y,Z +1), (X + l/2,r,Z), (X,r,Z), and 
{X,Y,Z) and t_hose of [TTTT] _arc {X ± 1/2, F 1/2, Z), 
(X, y, Z -f 1), (X, Y, Z), {X, Y,Z + 1/2), and (X, Y, Z). 
From these generators we identify the space groups as 
indicated in Fig. [181 



B. The star of N 



Again the analysis parallels that for RP214. Now the 
analog of Table |W1 is Tabic EB 



TABLE VI: As TableU but for the RP327 variables of the star 
of N = (1/2, 0, 1/2) (shown in Fig. |9j. T = (1/2, 1/2, 1/2), 
= (1,0,0), and Ty = (0, 1,0). All variables are odd under 
T = (0,0,1). 





7^4 


rrid 


■niz 


T 




Ty 


4>x,i 


0!/.l 


<Py.l 


4>x,i 


<t>x,2 


~4>x,i 


0a;, 1 


<t>y,l 


—4>x,i 


(t>x,l 


't>y,i 


0y,2 


<t>x,l 


-0y,i 


0a;, 2 


<Py,2 


<Py,2 


-'Px,2 


~<Px,l 


— 0a;, 2 


0a;, 2 


0I/.2 


4>x,2 


4>x,2 


-<t>V,2 




0y,2 


-0y,2 
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X = z = out 

FIG. 18: (Color online) As Fig. [6] but for the star of X for 
RP327 (with invariance under z — >■ z + 1). x, y, and z are the 
tetragonal axes and and X, Y , and Z are the conventional 
lattice vectors after distortion. The filled magenta circle is 
the tetragonal origin. The outer + or — sign gives the sign 
of the z-component of displacement of the upper layer of the 
bilayer and the inner + or — sign gives the sign of the z- 
component of displacement of the lower layer of the bilayer. 

= [1111] for a) and [1111] for b). The new origin is in the 
z = Q plane. The third new axis vector is [001] in terms of 
the original tetragonal coordinates. 



The free energy for the star of N is then 

2 

fc=i 

+ f [4>l,i + 4>li + €a + <i>\^ + ^2. + F4 (52) 

This free energy is the same as for RP214, so the struc- 
tures for RP327 will be related to those of RP214. In 
Fig. [12] we show the structures in which the 0's for the 



upper layer are identical to those of the single layer in 
Fig. [TUl Note that of all the generators of Fig. (TUJ 
only inversion (X, F, Z) connects the z = n layers to 
the z = n + 1/2 layers. Also note that according to 
Fig. [17] inversion for n — 2 systems introduces an extra 
minus sign compared to the n = 1 case. This fact sug- 
gests that by appropriately placing the new origins (as 
we have done in Fig. [T^ one can pass from RP214 to 
RP327, otherwise keeping the structure (and the space 
group) unchanged. We thereby determine the generators 
to be the same as for the analogous structures in Fig. |TU] 
for \xyxy] the generators are {X + 1,Y, Z), {X, Y + 1,Z), 
{X, Y,Z + l),and_(X,Y,Z), which is Cj (PT=#2). For 
both [OlOT] and [llTT], the generators are {X,Y^Z_+1), 
{X -f 1/2, Y_+ 1/2, Z), {X + 1/2, Y - 1/2, Z), (X, F, Z), 
and {X,Y,Z) which is C^^ (C2/m = #12). Although 
these two structures belong to the same space group, 
they are different because their unit cells are different 
(see Fig. [TU]). 



C. The star of P 

1. 6 structures 

The situation for 6'-structures at the star of P is sim- 
ilar to that at the star of X except that the displace- 
ments change sign under z — s> z + 1. The structures 
with On, a = On,b havc the same symmetry as that shown 
in Fig. [TId: however, each single layer is replaced by 
a bilayer in which 9n,a = dn,b- So the space group is 
again I4i/acd=#142 (D^^). The structure with 0„^a = 
—9n,b and which changes sign under z -> z + 1 is il- 
lustrated in Fig. [T6b . This structure is generated by 
{X - 1/2, F-H 1/2, Z+ 1/2), (X_+l/2, r- 1/2, Z+ 1/2), 
{X + 1/2, Y + 1/2, Z - 1/2), {X, Y, Z), {X, Y,Z + 1/2), 
and (y + l/4,X + 3/4,Z+l/4) and is therefore I4i/acd 
~ #142 (D^^) coming from irrep P2. In identifying the 
generators we used Fig. [T71 



2. 4> structures 

To identify the space groups of the structures shown 
in Fig. [20] recall the discussion just below Eq. (|52]) . By 
shifting the new origin, we identify the generators as: 

(x,y,z),H (x,F,z + i/2), (x,F,z), (x-i/2,r + 

1/2, Z+ 1/2), {X + l/2,y ^1/2,Z_+ 1/2), and {X_+ 
l/2,r+l/2,Z-l/2)for a), (X-hl/4, r+1/4, Z), {X,Y+ 
1/4, Z-Hl/4), (X, r, Z), (X-1/2, y+1/2, Z+1/2), (X + 
1/2, r - 1/2, Z_+ 1/2), and_(X 1/2, y 1/2, Z - 1/2) for 

b),and (x,y,z), (x,y,z), (x + i,y,z), (x,y + i,z), 

and {X, y, Z -|- 1) for c). These lead to the space groups 
listed in the figure. Because of the Lifshitz instability the 
wave vector is an incommensurate one close to P, as is 
discussed in Sec. III. 2. 
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FIG. 19: (Color online) As Fig. [T8]but for the star of N (with 
sign change under z ^ z + 1) for a) \xyxy], b) [0101], and c) 
[TlTT]. The new origins are at z = 1/4 for a) and at z = 1/2 
for b) and c). 



FIG. 20: (Color online) As Fig. [T8]but foi^the star of P (with 
sign change under z ^ z + 1) for a) [1111], b) [0101], and c) 
[xyxy]. The new origins are at z = 3/4 in a) and in b) and at 
z = 1/4 in c). The new out-of-plane axes are (0, 0, 2)t. 



V. RP SYSTEMS WITH n > 2 

Now we are in a position to analyze the situation of RP 
systems with n > 2. For the 0-structures the important 
issue is whether the eigenvalue X{Mz) which gives the 
symmetry of the stacking sequence is +1 or —1. For the 
^-structures the important issue is, as stated in Ref. |20| . 
whether n, the number of layers per substructure, is even 
or odd. 



A. ^-structures 

We first consider ^-structures. The 6'-structures asso- 
ciated with the star of either X or of P arc governed by 
the free energy 

2 n 

m— 1 k—1 

2 n 

+ Ak,iem,k9m,l + V , (53) 

m— 1 k,l — l 

where V contains interaction terms (of order A'^) between 
the two different slabs. Note that quadratic terms like 
di,aS2.p are excluded because they are not invariant un- 
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Top view of 
layers #1 and #k 




FIG. 21: (Color online) The ^-structures for one of the n-layer 
slabs. Left: The cross section of the n-layer slab showing the 
{x, y, z) coordinate system used to discuss the symmetry of 
the 6'-structures. Right: Top view showing that the structure 
(red dashed cross section) is characterized by giving the value 
of the rotation 6 for each of the n layers. 



der 7?.4, as one sees from TablelVlor Fig. [TH The ordering 
vector of the first n-layer slab in the unit cell is 

*1 - [014,^^1.2,^1.3 •■•^l,n] • (54) 

Note that $ is proportional to the eigenvector of the ma- 
trix A which has the minimal eigenvalue (so that it is 
the one which first becomes critical as the temperature 
is lowered). The matrix A is invariant under the mirror 
operation which interchanges layers k and n + 1 — fc, 
so that Ok O On+i-k- Therefore we know only that the 
eigenvector is either even or odd {i. e. the eigenvalue of 
Mz, \{Mz), is either +1 or —1, respectively), depending 
on the details of the interactions in the system. For in- 
stance, for n = 1 $1 = [1] and is even under M^. For 
n = 2 the critical eigenvector is either [11] (which is even 
under Mz) or [iT] (which is odd under Mz). For n = 3 
the critical eigenvector is either [lOT] (which is odd under 
Mz) or [a/3a], where a and /? depend on the interactions 
within the 3-layer subsystem and this eigenvector is even 
under Mz- On the basis of symmetry we can definitely 
not posit any specific form for the eigenvector (for n > 2), 
as is done in Table XIII of Ref. !20l. In all these exam- 
ples the ordering vector for the second n-layer subsystem 
obeys $2 = where the indeterminancy in sign re- 

flects the by now familiar frustration of 6 structures. 

For n > 2 the state we call I, given by the critical 
eigenvector of A, will develop at the structural phase 
transition at a critical temperature we denote T/. For 
this state 0^ will not be independent of k unless there is 
an unusual accidental degeneracy and therefore its free 
energy will be of the form 

Fi - ]^{T-Ti)Q] + a\Q\ . (55) 

The quartic term must be of order A because &j, is not 
independent of k. This state competes with state II which 



satisfies (for all k and I) 

26f.i = ^£xx = ^^yy ■ (56) 

State I becomes critical at a temperature Tj which is 
higher than that, T//, at which state II becomes critical. 
The free energy of state II can be written as 

Fii = liT-Tii)Qjj + /3Qjj . (57) 

Minimization with respect to the order parameters yields 

1 (T ~ T )^ 

Fi = --^ r-^ T<Ti, Fi=0, T>Ti 

2 aX 

Fii = --^ T<Tii, Fii=0, T> 

and these are plotted versus T in Fig. [211 One sees that 
for n > 2 we have a more complicated phase diagram 
than for n = 1 or 2. When A is large, for a small range of 
temperature phase I is stable, but at lower T we arrive 
at phase II. For n > 3 it is possible to have more than 
one phase transition before ultimately reaching phase II. 

F 
A 



Tll Tj 




FIG. 22: (Color online) The free energies of phase I (dashed 
line) and phase II (full line), as given by Eq. (|58p . These 
curves are drawn for qA//3 = 9. 

The possible 6'-structures of phase I depend on the 
wave vector (X or P) and whether the interactions within 
the system select X{Mz) = +1 or —1. First consider the 
star of X, for which 9m+2,k = Om,k- If K^^z) = +1 i^^ 
in Fig. [7^), then, irrespective of the value of the number 
of layers per subsystem, we will have a structure similar 
to that of Fig. [7^ (in which each layer of Fig. [7] is re- 
placed by an n-layer slab) with space group Cmca=#64. 
If interactions select, A(Mz) — —1 (this is only possi- 
ble for n > 1), then we have a structure of space group 
Ccca=#68, similar to that shown in Fig. [T5k . Next 
consider the star of P, for which 9m+2,k = —Qm.k- If 
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(a) 



n = 8 M^ = +l 
• • • • _ 



n = 8 M^=-l 
• • _ 



(b) 



n = 7 M^ = +l 



• • • 



(c) 



-► k 



n = 7 Mj^ = +1 



(d) 



TABLE VII: Summary of results for commensurate structures 
for RP n layer systems. Q denotes the wave vector, Var la- 
bels the angular variable, and A is the eigenvalue of the mirror 
operation within the n-Iayer substructure, as discussed below 
Eq. (|54p . In the last column we indicate whether the tran- 
sition is allowed to be continuous (Y) or not (N) or whether 
there is a Lifshitz instability (L). See Ref. [l^. The results in 
this table can be compared to the results for the X point given 
in Table XIII of Ref. 20,. Our analysis only allows structures 
#2, 5, 9, and 12 of that reference. 



Var A 



X 
X 



-1 
+1 



+1 

-1 



X 
X 



-1 
+1 



N (/) -1 



N 



+1 



+1 



Space group(s) 



#68 (Dl) 



See Fig. Y,N,L 



Ccca n > 1* 



#64 {DID Cmca n > 1 



14a 
7a 



#142 (D^'^) I4i/acd n > 1 
#142 {-Dfl) I4i/acd n > 1* 



7b 
14b 



#64 (D^^) Cmca 
#66 (Dll) Cccm 
#63 (D^^) Cmcm 
#67 (Dli) Cmma 



odd 
odd 
even 
even 



8b 
16c 
16a 



#12 (C^O C2/m 
#12 (CiJ C2/m 

#2 id) PI 
#12 (CL) C2/m 
#12 (CL) C2/m 

#2 {C}) PI 



odd 
odd 
odd 
even 
even 
even 



10b 
10c 
10a 
17b 
17c 
17a 



#72 (D^rj 
#70 (Di^) 
#15 (CiJ 
#72 (Di«) 
#70 (Di^) 
#15 (CiJ 



Ibam 
Fddd 
C2/c 
Ibam 
Fddd 
C2/c 



odd 
odd 
odd 
even 



11a 
lib 
11c 
18a 
18b 
18c 



Y 
Y 



Y 
Y 



Y 
Y 
Y 
Y 



N 
N 
N 
N 
N 
N 



L 
L 
L 
L 
L 
L 



* If the dominant interlayer interactions are those between 
adjacent layers, then, as discussed in the text, n must be 
even for this case to occur. 



FIG. 23: Typical variation of 6 as a function of the layer 
index, k, for dominant nearest neighbor interlayer interaction, 
Ak,k+i = w. For case (a), n is even and w; < 0, so that 
X{Mz) = +1. For case (b), n is even and w > 0, so that 
X{Mz) — —1. For cases (c) and (d) n is odd. For case (c) 
w < and for case (d) w > 0. In both cases A(Mz) = +1. 
To obtain X{Alz) = — 1 for n odd requires further neighbor 
interlayer interactions. 



A(Mz) = +1 (as in Fig. [Tb), then, irrespective of the 
value of the number of layers per subsystem, we will have 
a structure similar to that of Fig. [7)d with space group 
I4i/acd=#142. If interactions select X{Mz) = -1 (this 
is only possible for n > 1), then we have a structure of 
space group I4i/acd=#142, similar to that shown in Fig. 

The above remarks relied only on symmetry. However, 
now we consider the likely form of the interaction matrix 
Ak^i in Eq. ([55)) which determines 6^ as a function of the 
layer index k. If the dominant intraslab interactions are 



those between adjacent layers, then, as illustrated in Fig. 
[231 we obtain configurations analogous to ferromagnetic 
(panel a) or antiferromagnetic (panel b) spin structures. 
Thus, with nearest neighbor interlayer interactions, if n, 
the number of layers per slab is even, we can have ei- 
ther sign of A(Mz). If n is odd, then this special ansatz 
of nearest neighbor interlayer interactions can only give 
X{Mz) = +1. If, experimentally, the case X{Mz) = —1 is 
observed for odd n, one could conclude the existence of 
significant longer ranger interlayer interactions. (Such a 
situation is obviously possible in the presence of Coulomb 
interactions.) 

Our results are summarized in Table [VTIl Note that for 
the 6 structures the controlling variable is not the number 
of layers n, but A(Afz). It is interesting to note that for N 
and P, the (/)-structures for even and odd X(Mz) are very 
similar. Apart from the fact that their substructures are 
different, they only differ in the location of the center of 
inversion symmetry. 
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VI. DISCUSSION AND CONCLUSION 

We did not deal with the positions of the ions at the 
center of the octahedra or those between the layers of 
octahedra. Here we are only concerned with such ionic 
positions as they are modified by the orientational struc- 
tural transition. Each such ion sits in a stable potential 
well. The question is whether or not for systems with- 
out any accidental degeneracy there is a bifurcation so 
that additional space groups could be allowed when the 
positions of these "inessential" ions are taken into ac- 
count. The stable potential well can be distorted and 
the placement of its minimum will be modified by the 
octahedral reorientation. But a single minimum of a sta- 
ble potential well can not be continuously deformed into 
a double well without assuming an accidental vanishing 
of the fourth order term in the local potential. Similar 
arguments show that the perturbative effect of the cen- 
ter of mass coordinates of the nearly rigid octahedra do 
not produce anomalous effects. Of course, parameters of 
Wykoff orbits which are not fixed by symmetry will be 
perturbatively modified at the structural phase transi- 
tion. Similarly, the elastic strain tensor will be pertur- 
batively modified consistent with the symmetry of the 
resulting phase at the transition. 

Experimentally, it is striking that the structures ob- 
served as distortions from the tetragonal phase are in 
our much shorter list. For instance, in the data cited on 
p 313 and ff of Ref. [13 five systems with tilts are shown 
which go into either Cmca (64) or P42/ncm (138), except 
for Rb2CdCl4 whose structure is uncertain: either Cmca 
or Fccm (which is on neither our list nor that of Ref. 
[TtI because it involves two irreps). Systems (other than 
Rb2CdCl4 subsequently discussed in Ref. [20) in Table 
III of Ref. Tt' likewise go into either Cmca or P42/ncm. 

To summarize: we have analyzed the possible struc- 
tural transitions of the so-called Ruddlesden- Popper per- 
ovskite structure (such as K2MnF4 or Ca3Mn207, etc.) 
using a variant of Landau theory in which the constraint 
of rigid oxygen octahedra is implemented and our results 
are compared to the well-known results of Refs. [l3 and 
[20I . A check on the accuracy of our treatment of sym- 
metry is that our list of allowed structures (for K2MgF4) 
which can be reached via a single structural phase tran- 
sition is a subset of the list of Ref. [l3- We find that the 
rigid octahedral constraint eliminates all the structures 
in Table I of Ref. [13 for which the octahedral tilting 
transitions are discontinuous. It is also appealing that 
structures which are allowed by symmetry but which in- 
volve undistorted sublattices are eliminated by the octa- 
hedral constraint. The results for the K2MgF4 structure 
are summarized in Tables ITll and ITlll where one sees that 
our list of possible structures is a much reduced subset of 
Ref. [l3- For these systems, our analysis allows (see Ta- 
bles |TT] and |in| only 13 of the 41 structures listed in Ref. 
[TtI and of these only nine are commensurate structures. 
A summary of our results for commensurate structures 
for A„+iB„C3„+i is given in Table [VlTj 
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Appendix A: Minimization to order A ^ 

Here we show that corrections of order 1/ A to the mini- 
mization of the free energy do not qualitatively affect our 
conclusions. We consider the free energy 



F 



cga A 



-01 

2 1 



^yy 



+ \b [t.^^ + tyy) ^- cc^,] [Bl + B, 



(Al) 



where all the coefficients except A are of order unity and 
the summation is in Voigt notation. We assume that 
nonlinear elastic terms (of higher than quadratic order) 
in e can be neglected. We now set 



1 



[Bl 



(A2) 



The terms involving A become 



F(A) = ^cga'X {Bl -Blf + 2cea'X[e + C,'] • (A3) 

Therefore the quartic terms in Bk are of the form of Eq. 
(j^ where v is surely negative for large A. Accordingly, 
to all orders in 1/A we may set 



9l = Bl EE 0iX) 



(A4) 



Also we see that the coefficients of terms linear in Cxx 
and Cyy are equal. So we may set = = S,- Then ^ is 
determined to leading order in 1/A by 



Acga^XS,^ - \-f\B^ + UB^ 



1 2 



where 



-[(46 - xii)e + 2ce,,]e^ + -C44ely 
~dexy9iB2 , 



U = i(w + w)-26+icii 



(A5) 



(A6) 
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Then, by minimizing F of Eq. (jA5l) with respect to 
we get ^ '-^ or 

Sxx = eyy ^ -^e^ + e\-^ , (A7) 



where e is of order A°. This tells us that the octahedral 
mismatch (or distortion) is not zero, but is of order A~^. 
The corrections of order 1/A to the elastic constants, do 
not, of course affect the symmetry of the structure. 
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